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C Abstract 

^ We develop a master equation formalism to describe the evolution 

of the average density matrix of a closed quantum system driven by 
a stochastic Hamiltonian. The average over random processes gen- 
^-H erally results in decoherence effects in closed system dynamics, in 

1^ addition to the usual unitary evolution. We then show that, for 

an important class of problems in which the Hamiltonian is propor- 
tional to a Gaussian random process, the 2nd-ordcr master equation 
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\^ yields exact dynamics. The general formalism is applied to study the 

examples of a two-level system, two atoms in a stochastic magnetic 
I field and the heating of a trapped ion. 

^ PACS number(s): 03.65.Ca, 03.65. Yz, 05.40.-a, 03.67.Lx 

• • 

1 Introduction 

Cd The density operator encapsulates all the statistical information about the state 

of a quantum system. The evolution of the density operator of a closed system 
is governed by the Hamiltonian. In practice, the Hamiltonian can seldom be 
strictly determined or precisely controlled it fluctuates both in a temporal 
sense and between repeated realizations, which can be mathematically described 
by random processes. Therefore, instead of treating any Hamiltonian as deter- 
ministic in an idealized manner, we would like to take such fluctuations into 
account explicitly when studying quantum dynamics. Our goal is to obtain the 
average dynamics in the following sense: Suppose an ensemble of systems are 
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prepared in some initial state and subsequently evolve under a randomly fluc- 
tuating Hamiltonian, how does the density matrix that describes the ensemble 
as a whole evolve? 

Previous work on stochastic average dynamics was done by Budini jT| using 
a variational calculus method and Novikov's theorem, and by Guha et al. [2] 
using a non-perturbative cluster cumulant method. A different kind of average 
dynamics over the time domain was studied by Gamel and James [3] , assuming 
the deterministic (i.e. non-stochastic) Hamiltonian but taking into account the 
finite time-window of measurements. The master equation formalism is also 
widely used in the study of open systems dynamics [1]. It should be noted that, 
despite the formal similarity, our study is on the dynamics of closed quantum 
systems and no environment is involved. 

In this paper, we will first adopt a series expansion approach and derive a 
time-local master equation that describes the ensemble-average dynamics of a 
general quantum system. The general formalism is then used to study a repre- 
sentative class of Hamiltonians obeying Gaussian statistics. Finally, we apply 
the master equation method to some physical examples and find interesting 
phenomena such as fiuctuation-induced decoherence and decoherence-induced 
disentanglement. Throughout, our results are compared to exact dynamics and 
the validity of the master equation approach is discussed. 

2 Theory 

2.1 Ensemble-average density matrix 

Consider a closed, but not isolated, system for which the Hamiltonian is de- 
termined by some classical stochastic quantity x{t). Suppose an experiment is 
carried out repeatedly with each realization labelled by fx. The evolution of the 
density matrix that describes the quantum system in the /i-th realization 

is governed by the Hamiltonian H^^{t) — H[x^^{t)], and is given by 

p^{t) = ij^{tM)p^u''\tM), (1) 

where po is the initial density matrix, which is assumed to be uncorrelated with 
x{t) and thus is the same in all realizations. The unitary evolution operator 
U^{t,to) obeys the equation of motion, 

th^^U^{t,to)^H^{t)U^{t,to). (2) 
The average density matrix p{t) is defined as follows, 

m^^^^j^T^p^it). (3) 

It can be shown that p{t) is Hermitian, positive and of unit trace, which is 
ensured by the properties of the individual density matrices (t) . Thus the op- 
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erator p{t) is indeed a physical density matrix, describing the average statistics 
of the ensemble of realizations as a whole. 

The equation of motion for p{t) is formally given by 

^ 1^/9 1 ^ 

(4) 

However, since the right hand side cannot be written as a function of p{t), the 
equation is not of a closed form and thus not very useful. With the goal of 
obtaining a closed equation for p{t), we resort to a series expansion approach. 

2.2 Series expansion of the evolution operator 

Following the standard recipe for perturbative expansion [5 , the unitary oper- 
ator U^{t,to) in a particular realization p can be written as a power series in A 
(a parameter controlling the "strength" of the Hamiltonian) : 

oo 

where 

K(t,to) = I, (6) 
t 

Uli{t,to) = ^jdt'H^{t')ij^_,{t',to),n^l. (7) 

to 

Thus 'p{t) can be expressed in terms of U{^{t,to) and A: 



= E^™^"(^'^") Po Y.^''UUt,to) = ^A^-£fc[po] = £[pol (8) 

\ m / \ n J k=0 

where ffc[po] is the time-dependent map defined as 

k 

£k[po] = Yl^k-jit,to)poU]{t,to), (9) 

and £[pa] is a completely positive linear map [B]. Although their argument is a 
density matrix in this instance, £k and £ can act on any operator in general. 

2.3 Inverse transformation 

The map £ is a linear transformation that maps po to p. Since both po and p 
are operators in the same Hilbert space, and thus of the same dimension, it is 
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natural to postulate that an inverse transformation T = E ^ exists that maps 
p to po- That is, 

p,, = e-^\p\^T\p\. (10) 

Note that the meaning of "inverse" is purely mathematical here: the map 
is not to be confused with an inverse dynamical evolution in the physical sense. 
According to the semigroup property, a completely positive, trace preserving 
(CPTP) linear map is physically invertible if and only if it is an unitary map 
(see Section 3.8 of Ref.jTj). Therefore, in general, a CPTP map E : po does 
not have a physical inverse, that is, we cannot find another CPTP map that 
gives p — po- However, the mathematical inverse E~^ that serves our purpose 
here needs not be CPTP. 

Since the composition of a transformation and its inverse is the identity 
transformation, the identity J^[£[p]] = I[p\ holds for an arbitrary operator p. 
Following [3 , we adopt the ansatz that T can be expanded in powers of A, 

oo 

Hp] = E A™J"„[p]. Then we have 

oo oo oo I k \ 

^ A"j:„,[^A"f„[p]] = Y.^^ Y.H[^^-M] = A"^[p]. (11) 

m=0 ri=0 fe=0 \J=0 J 

Collecting terms of like powers in A, we obtain the set of equations involving 
J- and En- 



M£o[p]] - i[p], (12) 

M^iip]] + H[Eo[p]]^0, (13) 

M£2[p]] + Ti[Ei[p]] + T2[Eo[p]] ^ 0, (14) 

and so on. Solving for J-m in terms of En, and making use of Eg = I as defined 
m Eq.(|9|, we have 

To[p] = E„[p]=I[p], (15) 

= -£i[p], (16) 

^2[p] = -E2[p]+E,[EMl (17) 

and so on. 

2.4 Master equation 

Differentiating Eq. ([s]) with respect to time and making use of the inverse relation 



in Eq.(lO), we obtain the following equation: 



^^l^^p{t) = *M[po] - ihE\T\-pm- (18) 
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Here, the notation £[p] means first taking time-derivative of the time-dependent 
transformation £ to obtain a new transformation denoted by £ and then let- 
ting £ act on p; the argument p is not differentiated whether or not it is 
time-dependent. Assuming the order of differentiation and summation can be 
switched, we have 



k 

£{mm = Y.^"£n[J2\'^j',n[pim = j^x' | Y.£ATk^Mm | ; (19) 

m=0 fc=0 \j=0 

thus the equation of motion can be written as 
dt 



- E^' I ^f^T.^A^'^-Mm I ^ Y.^'c.ipit)]- (20) 

fe=0 \ j=0 J k=0 



Evaluating J-m and £„ explicitly, we find 



£o[p] = ih£o[To[p]] = 0, (21) 
Clip] = ih£o[Ti[p]]+ih£i[To[p]] = Hp-pH, (22) 
C2[p] = ih£o[T2[p]] + ifi£i[J^i[p]] + ih£2[Mp]] 



HUip -HUip + HpUl - HpUl 



-pUlH + pUlH -UipH + UipH, (23) 

and so on. Note again that the argument p is not to be averaged or differentiated 

and that terms like H are time-dependent just as Ck[p] are time-dependent 
transformations. 



Keeping terms up to 2nd order and setting A = 1 in Eq.(20), a time- local 
master equation is thus obtained for the evolution of p{t): 

ih^^m = [S,pit)]+Mt)-pm^ +'D[p{t)], (24) 



where A = HUi - HUi and V[p] = HpUl - HpUl - UipH + UipH. The 
effective Hamiltonian responsible for unitary evolution is 

H,ff^H+\{A + A^), (25) 
with which the master equation can be written in a more insightful way. 



Pit) = [Heff,m] + \ [A-A\-p{t)}+'D[m]- (26) 

It can be shown that the right-hand side of the equation can be put into the 
Lindblad form, which ensures Hermiticity, complete positivity and trace preser- 
vation of the evolution. 
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This result is formally similar to a previous work on average dynamics [5]. 
However, the physical meaning is different since the derivation in that case is 
for a time-average density matrix in a single realization. Incidentally, our result 
may also be reminiscent of some master equations for the reduced density matrix 
of open quantum systems. But it should be emphasized that our derivation is 
for a closed system and thus quantum entanglement with environment does not 
play a role here. 

Note that the above results are formally applicable to an interaction-picture 
density matrix, though we implicitly assume the SchrA^dinger picture in the 
derivation. The only difference is in the interpretation of the density matrix: 
When we use the interaction-picture density matrix pj, the expectation value 

of an observable O is given by (O) = Tr ^O/p/^ , where Oj is the interaction- 
picture operator instead of the original operator in SchrA^dinger picture. 



3 General applications 

3.1 Time-independent Hamiltonian 

Let us first apply our general result to the simple case where the parameters 
in the Hamiltonian are time-independent. That is, H — H{a). where a repre- 
sents random variable(s) instead of random process(es). Suppose further that 
H is of zero-mean, which implies a particular choice of "picture": Any time- 
independent, deterministic part of the Hamiltonian plus the average compo- 
nent of the stochastic part can be removed by a gauge transformation, that 
is, by switching to an suitably chosen interaction picture [8J. Note that, in 
the case of time-independent random variables, C/i(t, io) = {t ^ to)H/ih and 
Ul{t,tQ) = -C/i(i,to), thus [Ui, H] = [Ul, H] = 0. So we have i it = q and 
thus H^ff = 0. This result is special to the time-independent case, however. As 
we will see later, the effective Hamiltonian (to 2nd order) is in general non-zero, 
due to the non-commutativity of Ui and H in the time-dependent case. After 
simplification, the 2nd-order master equation in this particular case is 

^-p{t) = -L{j2-p{t)} + ^Hpm. (27) 
A class of problems of physical interest has a Hamiltonian of the form 

H = hY^aJn + alhl (28) 

n 

where a„ are jointly circular complex Gaussian random variables of zero meanj^ 
Substituting Eq.(|28| into Eq.([27|, we fincQ 



^According to the central limit theorem, Gaussian statistics is applicable when the random 
variables are due to the addition of many uncorrelated random sources. 

■^For those readers who might be concerned about the factor t on the right-hand side of 
Eq.l|29f : It is just a result of a time integral, as can be seen in the more general case of 
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—pit) = tY^iVki [-huh\p{t) - p{t)hkh] + 2h]p{t)hk) 

k,l 

+ (-hlkm - p{t)hlk + 2hip{t)hl)}, (29) 

where Tki — aka*i are the correlation functions. Note that the equation is of 
the famihar Lindblad form, which can be further simpHfied to a diagonal form 
through a linear transformation of the coefficients. 

3.2 Single real Gaussian random process 

Now consider the case of a time-dependent Hamiltonian 

H{t) = ha{t)h, (30) 

where a{t) is a (real) Gaussian random process. This is representative of a wide 
class of problems, for example, the Zeeman effect, where a{t) is proportional 
to the external magnetic field and h is the ^-component of the total angular 
momentum [5]. The random process in this section is taken to be the most 
general case, that is, we do not assume any additional property like zero-mean 
or stationarity. 

The ensemble-average dynamics under this Hamiltonian is exactly solvable, 
so let us first work out the exact, analytic result. The unitary evolution operator 
in a particular realization p is 

f7^(i, to) = exp (^-iv^{t)hj , (31) 

where 

i 

w^(i) = J dt'a^'it'). (32) 

to 

For an initial state p(to) = ^Pki{tn)\k){l\, where {|n)} is the energy-eigenbasis 

k,l 

with h\n) = E„\n), the evolution in a particular realization is 

p^^{t) = U^\tM)p{h)U'-\tM) - 5I/'«(^o)exp{-*«^W(^fc - El)} \k){ll 

k,l 

(33) 

thus the ensemble-average is 

time-dependent Hamiltonians later. Note that the trace-preserving property of the equation 
is guaranteed by the Lindblad form. 
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m = ^Pfc;(to)cxp{-w(<)(£;fe - Ei)}\k){l\. (34) 

k,l 

Invoking the special properties of Gaussian statistics]^ it can be shown that 



e^^{-w{t){Ek - El)} = exp <^ -i{Ek - Ei)v{t) 



(Ek - Elf 



v{tY - v{t) 



(35) 

Thus, the exact ensemble-average dynamics is given by the elements of the 
average density matrix: 



(36) 



-Pkiit) = Pki{t,)e^T^{~i{E,^Ei)v{i^^^^^-^^[v{tY~v{tr]\.{il) 



Now let us solve the same problem by the master equation approach. Using 
the results from Eqs.( 25|26 l, the following expression is obtained, 



i;i-p(i) = [fm{t)h,p{t)] 



a{t) J dt'a{t') - / dt'a{t)a{t') 
to 



2ih 



(h^ Pit) + Pity 



J dt'a{t)a{f) - a{t) j dt'a{t') 

to ^0 



hp{t)h, (38) 



which can be simplified to 
d 



dt 



p{t) = -ia{t) [h, p{t)] + D{t) h, [h, p{t)] 



(39) 



where 



D{t) EE a{t) / dt'a{t') ~ / dt'a{t)a{t') 



(40) 



To find the solution to this differential equation, we first write it down in terms 
of the matrix elements in the /i-eigenbasis: 

■'Since v{t) is a linear filtered Gaussian random process, it is a Gaussian random process 
itself. (See page 83 of Ref.[l0].) 
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dt 
dt 



Pkkit) = 0, (41) 
p,i(t) = \-i^t)(Ek~Ei) + D{t){Ek-Eif]p^i{t), (k^l). (42) 



Now we have a set of (de-coupled) linear ordinary differential equations (ODE's), 
which is easily solvable, 



Pkk(t) = Pkkito), 



(43) 



Pkiit) = Pki{to)exp { J dt' [~iait'){Ek - El) + D{t'){Ek ~ EiY 

puito) exp {-i{Ek ~ Ei)W) - ^^"'^'^' [^ - \ -(44) 



The second equality in Eq.(44) is obtained after some calculation, where v{t) is 
given by Eq. ( 32 1 . Thus we find the dynamics generated by the 2nd-order master 



equation coincides with the exact dynamics in this case. 

We observe that the energy population is conserved during the evolution 
while the coherence between different energy levels decays. Thus the evolution 
of the average dynamics is pure decoherence, with the "pointer basis" [12^ being 
the energy-eigenbasis. Note that, although the Hamiltonian varies with time and 
across different realizations, the energy-eigenbasis is the same throughout. In 
the case where some energy level is degenerate, we readily have a "decoherence- 
free subspace", in which quantum information can be stably stored [13]. Inci- 
dentally, a derivation in the context of open systems also suggests that energy 
eigenstates can emerge as "pointer states" in the so-called "quantum limit of 
decoherence" |15;. In that case, however, the decoherence results from small en- 
vironmental perturbation to the system, not from the fluctuation of the system 
Hamiltonian itself. 

We could have worked out the higher-order terms (i.e. >Cn[p] for n ^ 3) 
explicitly to see how accurate the 2nd-order approximation is. However, since 
the solution to the 2nd-order master equation coincides with the exact dynamics, 
we can readily conclude that all higher-order terms must sum up to zero without 
actually carrying out further calculations. 



Note that when deriving Eq.( 39 1 we do not assume anything about the nature 



of the random process a(i), not even the Gaussian statistics. In other words, the 
solution to the 2nd-order master equation is given by Eqs.(|43]|44| in all cases. 
On the other hand, the exact dynamics Eqs.( 36]|37 1 is based on the assumption 
of Gaussian statistics. If a{t) is not a Gaussian random process, then the exact 
dynamics will be different The implication is that, for a(t) being non-Gaussian, 
the 2nd-order master equation is not exact. 



''it may not be exactly solvable, but we know for sure that the solution is different from 
that in the Gaussian case. 
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3.3 Multiple jointly circular complex Gaussian random 
processes 

Let us briefly present the results for the more general Hamiltonian H(t) = 
fi,^ ( a„(t)/i„ + a* (t)/ijj ] , where a„(t) are jomtly circular complex Gaussian 

random processes of zero mean. The 2nd-order master equation in Lindblad 
form is found to be 



k,l k,l 

X i^-hkhjpit) - p{t)hkh\ + 2h\p{t)hk - h\hkp{t) - p{t)h\hk + 2hkp{t)h\'^ , 

(45) 



where 



aki{t) = ^ J dt' (ak{t)a*{t') - al{t)ak{t')) , (46) 

to 

t 

Pki it) = l [ dt' (afc(t)ar(i') + a*(tK(F)) ■ (47) 



to 



By comparing with Eq.(29l for the time-independent Hamiltonian case, we no- 



tice a major difference in this case is that the effective Hamiltonian is non-zero 
despite H{t) = 0. This effective unitary evolution results from the fact that 
the Hamiltonian operators at different times do not commute with each other 
in general. 

The 2nd-order master equation yields exact dynamics only for the special 



case of a single real Gaussian random process. In this more general case, Eq.(45 
does not lead to exact dynamics in general. This can be shown by explicitly 
evaluating higher-order terms like C/^lp] to flnd that they do not vanish in gen- 
eral. Despite the lack of perfect agreement, the master equation is nevertheless 
of great use in such cases, because the exact dynamics is generally not obtain- 
able and the 2nd-order master equation serves as a good approximation when 
the higher-order terms (e.g. /34[p]) are small compared to £2[p]- 



4 Physical examples 

We will illustrate the general results by applying them to a few examples. The 
findings will then be used to gain physical insights, and the validity of the master 
equation approach will be examined by comparing to the exact dynamics. 
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4.1 Two-level system 

First consider a two-level system (e.g. spin-1/2) subject to the Hamiltonian 
H{t) — huj{t)Z, where Z is the ^-component of Pauli operator and uj{t) a 
stationary Gaussian random process of zero mean. This falls into the category 



of Hamiltonians iSOl. Using Eq.(39), the 2nd-order master equation is obtained, 



= Z,[Z,p{t)] , (48) 



where d{t) = 4 J^^ dt'uj(t)uj{t'). Assuming an auto-correlation function of the 
form uj{t)uj{t') = ojg exp (— |i — t'l/T), where luq = oj{t)^, we have d{t) = 
4^T (1 - e-(*-*«)/'^) for t > to- 

Written in the Z-eigenbasis {|0), |1)}, the master equation becomes a set of 
linear ODE's: 

^Poo(i) = 0, (49) 
^PiiW - 0, (50) 
^PoiW = -d{t)Poiit), (51) 
^^Pwit) = -d{t)Pwity, (52) 

the solutions to which are 

Poo(*) = Poo(^o), (53) 
Pii(i) - Pii(io), (54) 

Poi(i) = poi(io)exp|-4air2(^^;^+e-(*-*")/^-l)|, (55) 
-p,,it) p,o(i„)exp|-4^r^(^^+e-(*-*o)/T_i^| (5g) 

As already discussed in the general case of a single Gaussian random process, 
the energy population remains constant while the coherence decays. This can 
be understood from a more physical perspective. Quantum coherence depends 
on the relative phase between the two components |0) and In an individual 

realization, the relative phase factor is exp J^^ (ii'w(t')|. Since u!{t) is a ran- 
dom process, the quantity dt'uj{t') becomes increasingly randomized with the 
passage of time. When the average is taken over an ensemble, these randomly 
distributed relative phase factors cancel out, thus suppressing the coherence. 
This suggests that quantum interference is difficult to observe because random 
fluctuation is ubiquitous. 
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As has been shown in the more general case Eqs. 
ter equation gives exact dynamics. When we work out the exact dynamics 
directly, the result is indeed found to be consistent, though such a direct calcu- 
lation is more demanding. Clearly, calculational convenience is one advantage 
of the master equation approach. 



43p4 ) , the 2nd-order mas- 



4.2 A pair of two- level systems in magnetic field 

Next consider an example of two atoms in magnetic field, each atom being 
a two-level system. The interaction of the spin with the B-field is given by 

H{t) = Hn{t) [z^ ® /" + ® Z"^ = hn{t)Ztotai, where Z^ is the usual PauU 
z-operator of the j-th atom. Suppose that the frequency fl{t), which is propor- 
tional to the B-field strength, is a stationary Gaussian random process of zero 
mean and that its auto-correlation is of a Markovian type il{t)il{t') — ^6{t—t'), 



where the constant 7 has dimension of inverse-time. Applying Eq.(39l, it can 
be shown that the 2nd-order master equation is 

d _ If- - - 1 

^P(^) = -J^l yZtotaU [Z total, P{t)]\^ ■ (57) 

Suppose the system starts in an entangled state between two atoms |5'(io)l 



^ (|01) -I- |10)). Since it is an eigenstate of Ztotai, the right-hand side of Eq.(57l 
is identically zero, thus the system does not evolve (except possibly to an unob- 
servable global phase). So the two atoms remain entangled over time. Indeed, 
notice that |01) and |10) are degenerate eigenstates with the same energy. Thus 
any arbitrary superposition state of |01) and |10) will remain unchanged over 
time; in particular, the coherence between them does not decay. Thus, any state 
in this degenerate subspace is immune to decoherence, making it a good place 
to store quantum information |14| . 

Let us see what happens if the 2-atom system starts in a different entangled 
state like \'ii{to)) = ^ (|00) -I- |11)). Writing the master equation in the energy 
eigenbasis {|00), |01), |10), |11)}, we obtain a set of decoupled linear ODE's as 
usual. The solutions are found to be 

Poo,oo(i) = Piiai(0 = \: (58) 

PooaiW = Pii.oo(^) = ^exp{-7(t-to)}, (59) 

while the rest of the matrix elements are identically zero. Note that decoherence 
occurs here, as is expected, since the initial state does not lie in the decoherence- 
free subspace. Furthermore, as i — > 00, the coherence is suppressed to zero and 
m ^ i|00)(00| + i|ll)(ll| = i|0)(0| ^ |0)(0| + i|l)(l| ® |1)(1|. Interestingly, 
the two atoms become disentangled, as there is no quantum correlation between 
them. In contrast to the general belief that entanglement leads to decoherence, 
as is widely studied for open quantum systems, here we find that decoherence 
can result in disentanglement in the case of a closed system. 
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4.3 Heating of a trapped ion 



Consider an ion with mass M and charge e in a harmonic binding potential 
with characteristic frequency ojq. The ion is driven by a classical electric field 
E{t), which is a stationary Gaussian random process of zero mean. It is more 
convenient to work in the interaction picture, in which the easily solvable, de- 
terministic evolution induced by the harmonic potential is treated separately. 
The interaction-picture Hamiltoniarj^ is given by H{t) = ih [u{t)a^ — u*(t)a], 
where u{t) — ieE{t)e^'^°'^ / \/2Mhu!o and a (a^) being the zero-time annihilation 
(creation) operator for the harmonic oscillator. The evolution is exactly solvable 
and the analytic results are given in [IT,. 

Let us derive the 2nd-order master equation for this case. Note that it does 



not fall into the category of a single real Gaussian random process as in Eq. ( 30 1 



Since H = 0, applying Eq.(25l, the effective Hamiltonian is found to be 



t 

1 

to 

^2 



^ I" dt' (u{t)u* {t') - u* {t)u{t')j [a,a'f] 
t 

J dt'E{t)E{t') sin [ujo{t - t')] i. (60) 



e 

2Mwo 



In this case, the 2nd-order contribution to the effective Hamiltonian is non- 
zero, a consequence of the non-commutativity of H{t) and Ui{t). However, 
since H^f f is proportional to /, the unitary part of the equation of motion 
results only in an unobservable global phase in this case. For more general 
cases, Hef f can be different from the identity / and can well lead to non-trivial 



dynamics. Evaluating the remaining terms in Eq.(23l for this example, we find 
the following master equation: 



d 

—p(t) = ~C(t) (a)ap{t) + plt)a)a - 2ap{t)a}) 
at ^ ' 

- C{t) {da^p{t) + p{t)aa) - 2d^p{t)d) 

- e"*""* [C{t) ~ iSit)] {{a^fp(t)+p{t){a^f - 2a^p{t)a^) 

- e^2'"°* [C{t) + iSit)] {a^p{t) + p{t)d^ - 2dp{t)a) , (61) 

where C{t) {S{t)) are proportional to the incomplete cosine (sine) transform of 

^Throughout this section we work in the interaction picture. The subscripts to denote 
interaction-picture operators are dropped for notational simphcity. 
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the field correlation function, viz 

t 

e2 



^^^^ ^ ^^^^ / di'i?(i)i?(i')cos[^o(t-i')], (62) 



to 



'^^^^ ^ 2mu7^ I ^^'^(*)^(^') N(i - t')] . (63) 



to 



Assuming E(t)E{t') = £;(0)2 exp - t'|/r) and setting to = for conve- 
nience, we find C(t) = (l/2ri) {e~*/-^ [a;oTsin(wot) — cos(a;ot)] + l} and Sit) = 
— (1/2ti) {e~*/-^ [sin(a;ot) + a;oT'cos(a;ot)] — i^oT}, where ti is the heating time 
defined as l/n = (e^iJ^Op/Af/iwo) (T/(l + w^T^)). 

Unlike the case of Eq.(30l, the 2nd-order master equation does not generate 
exact dynamics in this case. To get an approximation of the heating from the 
ground state (i.e. pqo{0) — 1) for a short period of time, let us write the master 
equation in the energy eigenbasis of the harmonic oscillator. 



^Poo(t) - -2C(t)poo(t) + 2Cit)p,,it) 

- V2e2^""* [C{t) - iS{t)] po2(t) - V2e-2-ot [C{t) + iS{t)] p^^it)- (64) 

Since Puit), Po2(0 P2o(0 ^^e all negligibly small for t 'ti T, 1/wo, we 
have ^Poo(^) — ~2C(t)poo('^) to the lowest order. In the same manner, since 
the depopulation 1 — Poo(0 perturbatively small, we have Pog(t) = 1 to the 
lowest order on the right-hand side. Thus an approximate differential equation 
is obtained as ^Poo(0 — — 2C(t). Solving this ODE, we find, to lowest order|^ 



l-PooW-2/.t'C(t')^^t^ (65) 



which holds for short times and agrees with the analytic result in 

To investigate the evolution of the system for longer times, we write down 
the master equation in the same basis and solve it numerically. Since the Hilbert 
space is of infinite dimensions, it is not possible to write down the complete set 
of ode's for the matrix elements. Instead, we truncate it to a set of 5 x 5 
coupled ode's that includes only the matrix elements of the five lowest energy- 
eigenstates and their coherence The numerical solutions of F{t) = Poo(^) 

^The same result can be obtained by working out the evolution of Pii(t) for t <^ T, I/ujq 
using approximation to the same order. 

'^Since the system is of continuous nature, we could have solved the master equation in 
the Wigner representation. However, that approach is not analytically solvable either and is 
not computationally economical. Furthermore, even in that case, we still have to accept the 
imperfection of truncation since the numerics can only be done on a finite region of the "phase 
space". 
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Figure 1: The fidelity of the ground state as a function of dimensionless time ujot- Dash lines 
represent our numerical results, while solid lines are exact dynamics from 



(i.e. fidelity of the ground state) are shown in Figure 1 for different sets of 
parameters. It can be seen that, as ujqTi (i.e. the dimensionless heating time) 
increases with ujqT (i.e. the dimensionless coherence time of E{t)E{t')) fixed, 
the numerical result gives better approximation to the exact dynamics. Also 
note that, for larger values of wpT, the ground state population shows temporary 
revival against its general trend of decrease. 

Despite the artificial defect caused by the truncation of the set of ODE's, it 
is of more interest to know the validity of the 2nd-order master equation itself 
in approximating the exact dynamics. This is done by comparing the size of 
higher-order terms to that of the 2nd-order term. Using the Gaussian moment 
theorem 1 10|, it is easy to show that >Cn[p] = for all odd numbers n, so we are 
interested in the ratios between the even-number-order terms. Assuming loqT 
is fixed, it can be shown that £4[p] oc l/rj^wo as opposed to C2[p\ oc 1/ti, so 
Ci[p\/ C2[p\ oc I/wqTi. The same ratio holds for £6[/o]/'C4[/9], etc. Therefore, 
as long as I/wqTi is small, the higher-order terms become progressively small, 
lending legitimacy to the 2nd-order approximation. This is also consistent with 
the previous observation from the numerical results. Physically, this can be 
better understood by switching to the SchrA^dinger picture: The external field 
H field oc 1/ri is treated as a perturbation to the self-Hamiltonian of the sys- 
tem Hseif oc loq. Naturally, as the relative size of the perturbing Hamiltonian 
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Hseif /Hfieid oc 1/ujqTi becomes smaller, a perturbative method such as the 
2nd-order master equation gives better approximation to the exact dynamics. 



5 Conclusion 

In this paper we have presented the derivation of a master equation for closed 
systems driven by stochastic Hamiltonians from an ensemble-average perspec- 
tive. The principal result is given in Eqs.( 25][26 ). The validity of this approach is 



examined and 2nd-order master equation is found to yield either exact dynamics 
or good approximations to exact dynamics. 

Applying the formalism to various physical examples, we find the ensemble- 
average dynamics usually contains decoherence terms in addition to the unitary 
evolution. Decoherence plays an important role in the foundational problems 
of quantum mechanics, as it gives insights in two aspects of the measurement 
problem, namely the absence of observable superposition and the problem of 
preferred basis |12| . Extensive research has been done on how environmental 
entanglement causes decoherence in open systems. However, as our findings sug- 
gest, decoherence could also be attributed to the random fiuctuations of physical 
quantities in closed systems. If this is true, then the tension between the classi- 
cality of our experience and the quantumness of the underlying laws of physics 
could be reconciled in some degree by the ubiquitous random fiuctuations. Fur- 
ther investigation is needed to find out (a) to what extent decoherence is actually 
caused by random fiuctuations and (b) whether/how we can distinguish it from 
the usual entanglement-induced decoherence through physical observation. 
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